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We develop an adiabatic formalism to study the Hawking phenomenon from the perspective of 
Unruh-DeWitt detectors moving along non-stationary, non-asymptotic trajectories. When applied 
to geodesic trajectories, this formalism yields surprising results: (i) though they have zero accel- 
eration, the temperature measured by detectors on circular orbits is higher than that measured 
by static detectors at the same distance from the hole, and diverges on the photon sphere, (ii) 
in the near-horizon region, both outgoing and incoming modes excite infalling detectors, and, for 
highly bound trajectories (E <C 1), the latter actually dominate the former, (Hi) in this region, the 
relationship between the temperature of Hawking radiation and the relative velocity between the de- 
tector and the hole is not of Doppler type. We confirm the apparent perception of high-temperature 
ingoing Hawking radiation by infalling observers with E <C 1 by a flux computation. We close by a 
discussion of the role played by spacetime curvature on the near-horizon Hawking radiation. 
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I. INTRODUCTION 

Hawking has famously predicted 0, [2J that the grav- 
itational collapse of a (say spherically symmetric) mat- 
ter distribution of mass M will be perceived by ob- 
servers at future infinity in the form of a stationary, out- 
going, thermal flux of massless particles with tempera- 
ture Tjj = (87rM) _1 . Like most general- relativistic ef- 
fects, this "Hawking phenomenon" is a priori observer- 
dependent, and it is not immediately clear what other, 
non-asymptotic observers would measure. To investigate 
this question, Unruh introduced in [3j a most useful and 
practical approach, based on a simple particle detector 
model (referred to as an "Unruh-Dewitt (UDW) detec- 
tor"). In the case of static UDW detectors at "radius" 
r in a Schwarzschild spacetime, for instance, one finds 
[4| that the temperature measured by UDW detectors 
isT#(l — 2M/r)~ 1 / 2 ; this shows that Hawking radiation 
satisfies the Tolman equilibrium condition [5(, just like a 
normal thermal bath. 

Another interesting class of trajectories worth investi- 
gating are infalling geodesies. Since these are not tangent 
to a timelike Killing field, however, one cannot use any 
Tolman-like a priori argument to infer the temperature 
measured by UDW detectors along these trajectories, and 
in particular when they cross the Schwarzschild horizon. 
Would they record Hawking radiation there? Or would 
they not, because Hawking radiation is created at some 
distance away from the black hole [a]? 
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Unruh gave an answer to this question in [3]. Tak- 
ing his cues from the geometric similarity between the 
Schwarzschild and Rindler horizons, and the observation 
that the Unruh temperature Tjj = ajl-K of geodesic ob- 
servers in Rindler spacetime is zero (because their ac- 
celeration a is zero), he argued that "a geodesic detec- 
tor near the horizon will not see the Hawking flux of 
particles" . This conclusion — although not supported by 
an explicit calculation in [3j — fits with the general view 
put forward in that paper that, near the horizon, the 
in-vacuum of gravitational collapse is not different from 
the Minkowski vacuum. The reduction of the (curved 
spacetime) Hawking effect to the (flat spacetime) Un- 
ruh effect in the near-horizon region [7j is further sup- 
ported by the following fact: the Hawking temperature 
Tjy(l — 2Af/r)~ 1 / 2 perceived by a static observer at a ra- 
dius r > 2M in Schwarzschild spacetime approaches the 
Unruh temperature (Af/r 2 )(l - 2Af/r) _1 / 2 /27r for a ob- 
server with the same acceleration in Minkowski spacetime 
when R — > 2Af. Singleton and Wilburn have described 
this fact by saying that "the equivalence principle is re- 
stored at the horizon" [J]. 

For all that, heuristic arguments as well as actual com- 
putations have recently challenged this conclusion. In the 
former category, Heifer has argued [8] that "the vicinity 
of a black hole is a region in which essentially quantum- 
gravitational, Planck-scale, physics must dominate"; for 
Almheiri et al. [9(, the principle of information con- 
servation and other quantum-information-theoretic con- 
straints suggest the presence of a "firewall" at the hori- 
zon. While interesting in themselves, these arguments 
are unfortunately too sketchy to lead to a definite pre- 
diction concerning the fate of horizon-crossing geodesic 



detectors. Following a more conventional semi-classical 
approach, Barbado et al. [10( have recently devised a 
framework, based on the notion of "effective tempera- 
ture" , to study the nature of Hawking radiation near the 
horizon more explicitly. (The "effective temperature" in- 
troduced there is meant to give an estimation of the re- 
sponse of UDW detectors at any given point along their 
trajectory.) Somewhat surprisingly, they find that the ef- 
fective temperature of detector dropped with zero initial 
velocity from infinity will rise, until it reaches the value 
ATh when crossing the horizon. This result appears to 
contradict Unruh's conclusion [3|, and does not accredit 
the notion that the Hawking effect reduces to the Unruh 
effect in the near-horizon region. 

Very appealing for its conceptual simplicity, the "effec- 
tive temperature" approach of [HJ has several technical 
limitations which could legitimate some skepticism about 
these conclusions. First, the actual response function of a 
trans-horizon geodesic detector is not explicitly evaluated 
in [n| , nor is any solid argument given to the effect that 
the "effective temperature" defined by these authors is 
indeed the one measured by an Unruh-DeWitt detector. 
Second, the framework used in [ICf, being entirely based 
on the retarded Eddington-Finkelstein time coordinate 
u, does not apply to the interior region of the black hole; 
hence one could speculate that the effect found in pj| is 
perhaps an artifact of the horizon being singular with re- 
spect to the u-coordinate. Third, the vacuum state used 
in the definition of the "effective temperature" function in 
|10J ] is not the standard Unruh vacuum, but a new, uncon- 
ventional state in the extended Schwarzschild geometry. 
One may fear that this new state is perhaps unsuitable to 
describe the physics of Hawking radiation. Fourth — and 
more crucially — the fact that a radially infalling geodesic 
detector clicks at the horizon, and even records a finite 
temperature there, does not mean by itself that it per- 
ceives "Hawking radiation". The local spacetime geom- 
etry is not stationary along a radial geodesic, and actu- 
ally changes more and more rapidly as the horizon (and 
then the singularity) is approached; it could be that the 
temperature measured at the horizon is due to this (triv- 
ial) curvature effect, and not to the peculiar structure of 
the vacuum which is the true origin of Hawking radia- 
tion. Rather than "does a near-horizon detector record 
a non-zero temperature" , the question which should be 
addressed is therefore "does a near-horizon geodesic de- 
tector record a non-zero temperature that cannot be ex- 
plained away by a time- varying gravitational potential" ? 

The purpose of this paper is to shed more light on 
these questions by reconsidering the response of geodesic 
UDW detectors during and after gravitational collapse. 
We address the four aforementioned concerns as follows: 

1. We define and study a suitably defined time- 
dependent "temperature" function, as in [l0| . but 



also related it explicitly to UDW response functions 
by an adiabatic expansion. 

2. We only use globally defined coordinates, regular 
on the horizon as well as everywhere else (except at 
the singularity). 

3. We avoid working with the extended Schwarzschild 
spacetime, where the effect of gravitational col- 
lapse must be mimicked by a suitable choice 
of vacuum state (e.g. the "Unruh vacuum" of 
Q or the "collapse vacuum" of [10]). Instead, 
we consider gravitational collapse geometries and 
the corresponding — uniquely defined — in-vacuum 
state. 

4. In a final section, we also consider generalized (non- 
Schwarzschild) black hole geometries, where the de- 
generacy between the surface gravity and the cur- 
vature scale at the horizon is lifted. This allows us 
to disentangle curvature effects and Hawking radi- 
ation at the horizon. 

Another assumption made in 10] which we relax in 
this paper is that only outgoing modes couple to UDW 
detectors. In fact, we will see that there exists a class 
of infalling trajectories for which the response of UDW 
detectors is actually dominated by incoming modes. This 
is a somewhat surprising of our analysis, which, to our 
knowledge, was not anticipated in the literature. 

The plan of the paper is as follows. In sec. HH we in- 
troduce the model of gravitational collapse used int he 
paper (a Vaidya ingoing shell), and discuss our "quasi- 
temperature formalism" to study the response of UDW 
detectors on arbitrary trajectories in that spacetime. We 
then apply this formalism in sec. IHII to various geodesic 
trajectories: circular orbits, radially infalling trajectories, 
and inspiral orbits. Our findings are compared with a flux 
computation in sec. IIVI In the next section, we address 
in more details the role played by the local spacetime 
curvature on the response of near-horizon geodesic detec- 
tors, by means of an artificial black hole model where the 
curvature vanishes near the horizon. Sec. IVII contains a 
discussion of our results and our conclusion. Details on 
adiabatic expansions of UDW-like response functions are 
given in the Appendix. 



II. THE HAWKING EFFECT ALONG 
GENERAL TRAJECTORIES 

The Hawking effect is the perception of an outgoing 
thermal flux at temperature Tjj = (87rAf) _1 by asymp- 
totic inertial observers at rest relative to a Schwarzschild 
black hole. In this section we introduce a formalism 
to study response of particles detectors moving along 
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FIG. 1. Definition of the (giobaily defined) "eikonal coor- 
dinates" (v+ , V- ) for sphericaliy-symnietric gravitational col- 
lapse. 



more general trajectories. This allows us to consider two 
aspects of the Hawking phenomenon which arc not of- 
ten described in the literature: the onset of black hole 
evaporation after the horizon has formed, and the time- 
dependence of the spectra recorded by non-stationary de- 
tectors. 

Throughout this paper, we shall make the following 
assumptions and approximations: 

• We only consider minimally coupled massless scalar 
fields. 

• We neglect the contribution of non-spherically sym- 
metric (I y^ 0) field modes as well as all backscat- 
tering effects, so that the dynamics of the field is 
effectively two-dimensional. 

• We neglect any backreaction of Hawking radiation 
on the background spacetime. 



A. Collapse geometry and null coordinates 

Our model gravitational collapse in this paper is the 
Vaidya ingoing shell. (Another collapse geometry will be 
discussed in sec. [Vj) As is well known, this spacetime 
consists of two patches separated by a null thin-shell: a 
flat region inside the shell, and a Schwarzschild region 
outside the shell. Its metric is conveniently written in 
Eddington-Finkelstein advanced coordinates (v,r) as 



where r s = 2M is twice the mass of the shell (which is 
located at v = 0), Q(v) is the Heaviside function and 
dfl 2 = d9 2 + sin 2 9d(p 2 is the standard angular metric. 
The surface gravity of the black hole is k = (2r s ) . 

To analyze the Hawking phenomenon in such a col- 
lapse model, it useful to introduce globally defined null 
coordinates (v+,V—) for the (v,r) sector of spacetime. 
These are constructed as follows. For each event x, con- 
sider the two null rays (incoming and outgoing) meet- 
ing at x, and let v + (x) and V—(x) be their respective 
Eddington-Finkelstein advanced time at past null infin- 
ity, with v-(x) < v+(x) = v(x) (see Penrose diagram in 
Fig. [1]). 

The physical interpretation of the (v + , w_ ) coordinates 
is tied to the eikonal approximation for the propagation 
of massless fields (j>(x), according to which 
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4>(v+{x),r) - 4>(v-{x),r) 



(2) 



In this approximation, the field at a given point x is ex- 
pressed as the superposition of two spherical waves em- 
anating from past null infinity J~: a convergent one, 
arriving directly from v+(v,r) (first term above), and a 
divergent one, arriving from v—(v,r) after a reflection off 
the origin (second term above), as in Fig. [TJ (For a har- 
monic mode e~ zulv on J~ , the phases of these two waves 
at a given point x are u>v±(x)). For this reason, we call 
(v+,V— ) the eikonal coordinates of x. 

In the case of the Vaidya metric (jTJ) , it is possible to 
give the (v,r) i-» («+,«_) mapping in closed form. By 
construction, we have w + = v, and, when («, r) is inside 
the shell, V—(v, r) = v — 2r; for a point (v, r) outside the 
shell (v > 0), we must integrate the equation ds 2 = 
from (y, r) back to the point (v — 0, r) where it meets 
the shell, and then write V-(v,r) = —2r. This gives 



V-(v, r) = 



v — 2r 



-2r, 



1 + W(6e s -" V ) 



for v < 
for v > 0, 



(3) 



ds 2 



(l-—e(v)\dv 2 + 2dvdr + r 2 dn 2 , (1) 



where 8 = r/r s — 1 and W(z) denotes the Lam- 
bert VF-function, defined as the principal solution of 
W(z)e w ^ — z. (For the reader's convenience, the graph 
of W(z) plotted in Fig. H) 

Since W(z) ~zasz-)0, the identity (|3|) immediately 
gives the equation of the event horizon as u_ = — 2r s . 
But, more importantly, it also reveals the peculiar dis- 
ruption of phase fronts induced by gravitational collapse: 
from ([3]) we see that v^ is constant where Se 5- ™ is con- 
stant, and in particular u_ ~ — 2r s where 5e s ~ KV <C 1. 
Due to the peculiar analytic form of the function 8 h-> 
8e 5 ~ KV for each v > 0, the region where Se s ~ KV -C 1 is 
very sharply defined. Inside this region, which expands 
away from the horizon as v grows (until it eventually cov- 
ers the whole spacetime), the phase of a vacuum fluctua- 
tion emerging from the shell after bouncing off its center 
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FIG. 2. The Lambert W-function, defined on [e 1 ,cxa) by 
W{z)e w{z) = 2 and such that W(z) ~ 2 for 2^0. 



Notice that, in the Hawking region, the eikonal coor- 
dinate V- takes the same functional form as the Kruskal- 
Szekeres U coordinate for eternal black holes. Defining 
indeed, a la Kruskal-Szekeres, U(v,r) = -2r s e~ KU< - v ^ 
with u(v, r) = v — 2r — 2r s log S, from see from © and 
the asymptotic formula W(z) ^ z as z ~ > that 



V-(v,r) ~ U(v,r) — 2r s 



(4) 



when Je 5 "™ <C 1. Since eternal black holes are not real 
objects, the relation (Ql can be thought of as providing 
the physical interpretation of Kruskal's and Szekeres' U 
coordinate: it captures the structure of outgoing phase 
fronts in the Hawking region of a (non-eternal) black hole. 



Unruh-DeWitt detectors 



is almost exactly equal to — 2r s ; this is where the Hawk- 
ing phenomenon takes place. We shall call it the Hawk- 
ing region. The level curves of v_(v,r) are represented 
in Fig. [3] in pictorial terms, Fig. [3] shows how gravi- 
tational collapse "opens up", or perhaps "breaks open", 
the in-vacuum. 

I 



Unruh-DeWitt (UDW) detectors [1, [TTJ are point-like 
monopole detectors, which measure the Wightman func- 
tion G(x,y) of the field along a given trajectory. Specif- 
ically, the response function of a UDW detector moving 
along the trajectory 7(7-) is, to first order in perturbation 
theory, 



/oo />oo 

du X (u) / dsx(u-s)e- tns G(-f(u),~f{u~s)). 
-00 Jo 



(5) 



In this expression, Q is the energy gap between the two 
stationary states of the detector, G(^f{u),^{u — s)) is 
the pull-back of the Wightman function to the detec- 
tor's worldline, and x( u ) is a non-negative "switching" 
or "window" test function. The introduction of a switch- 
ing function in the definition of the UDW response func- 
tion ensures that 1Z(£1) is well-defined and consistent with 
causality: choosing a switching function such that xi u ) — 
for u > t and u < t — At ensures that TZ(Cl) only de- 
pends on the time interval At in the past of r. Motivated 
by this causality issue, Svaiter and Svaiter introduced in 
12[ the simple ansatz xi u ) = B(m — r — Ar)8(r — u), 
which in the At — > 00 limit leads to a transition rate 
(r-derivative of the response function 1Z) 



£(n) = 2Re / dse- ins G(~f(T),j(T-s)). (6) 
Jo 

However appealing at first sight, this ansatz has two un- 
pleasant features, which are due to the fact that x( u ) = 
<d(u — r — At)0(t — u) is not a proper test function (it is 
not smooth): first, the standard — iO prescription for the 
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Wightman function yields non-Lorentz-invariant transi- 
tion rates [13j and must be replaced by a more compli- 
cated prescription |13l - [la |: second, the ultraviolet behav- 
ior (large Q limit) is qualitatively different on stationary 
and non-stationary trajectory [17|, LL8[ . Both these short- 
comings disappear if a smooth switching function is used 
instead. 



C. Wightman function in the (s-wave) in-vacuum 

As mentioned earlier, in this paper we shall restrict 
our attention on the s-wave sector (spherically symmet- 
ric field configurations) of the in-vacuum. In this approx- 
imation, standard arguments [19|, [20] show that G(x, y) 
takes the logarithmic form 

G{x,y) oc ]n((v + (x)-v + (y)-i0)(v-(x)-v-(y)-i0)y 

(7) 
Using the relation Q , this expression can be identified in 
the Hawking region with the one given in [3j to describe 
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FIG. 3. "Portrait of the vacuum": level curves of the eikonal 
coordinate V-(v,r) in the Vaidya spacetime. The dashed hor- 
izontal line is the collapsing shell, the dotted vertical line is 
the event horizon, and the region in the top-left corner is the 
"Hawking region" . Its boundary, defined by V- (v, r) = — 2r — e 
for some e, is very sharp, and may be thought of as a "wall". 



a radiating eternal black hole and known as the "Unruh 
vacuum" , 



Gu(x, y) ex In [{v(x)-v(y)-iO) (U(x)-U(y)-iO)j. (8) 

These states — the collapse vacuum and the Unruh 
vacuum (jHJ), in the s-wave sector — have a peculiar prop- 
erty, which is not shared by more general vacua: the in- 
coming and outgoing modes couple to UDW detectors 
independently. This means that the UDW response func- 
tion 1Z splits as 1Z = 1Z + + 1Z- , where 1Z± is given by ((5|) 
with G{x, y) replaced by 



In (v± (-y(u)) - v± (l{u)) - iO) 



(9) 



The thermal properties of TZ+ and 1Z- can thus be stud- 
ied independently, and are completely determined by the 
behavior of each of eikonal coordinates v± along the de- 
tector's trajectory. In particular, in the case of stationary 
trajectories leading leading to a thermal spectrum, the 
thermality (aka detailed balance) condition for R + and 
1Z- defines two temperatures T + and T_, by 



TZ±(-Q) 



- e n/T ± 



1Z±(Vl). 



(10) 



These temperatures need not coincide in general. This 
is clear in the asymptotic Hawking effect, where T + = 
and T_ = Th = k/2tt. 



D. Hawking temperature from the peeling of 
outgoing modes 

The mechanism responsible for the perception of Hawk- 
ing radiation has been described by many authors as an 
"exponential redshift" or "peeling" effect. After Q, this 
statement is usually expressed in terms of the canoni- 
cal mapping from J + to J~ which relates the retarded 
and advanced Eddington-Finkelstein times of a given null 
ray. From the perspective of UDW detectors, however, 
it is more natural to state the "peeling" condition in a 
trajectory-dependent way, as follows: if 7(1") is the tra- 
jectory of an asymptotic observer at rest relative to the 
hole, the quantity t>_(7(r)) satisfies 






(11) 



(From now on, we shall drop the explicit reference to 
7, the trajectory being understood from the context.). 
Integrating this condition gives 



. (t) -V-(t- s) oc e -«(^- s )/ 2 sinh(Ks/2). 



(12) 



If we plug this relation into © and perform the stan- 
dard residue integration [19j , we arrive at the celebrated 
thermal spectrum at temperature Th = k/2tt 



7t_(ft)oc 



1 



ft(e 27r °/ K - 1)' 



(13) 



(The ft — > divergence is of course an artifact of the 
two-dimensional approximation.) 

The main advantage of this trajectory-based formula- 
tion of the exponential redshift argument is that it ap- 
plies mutatis mutandis to any trajectory where u_ /t>_ is a 
constant: for the same reasons that asymptotic observers 
such that \v/v\ = k perceive a thermal spectrum at tem- 
perature k/2tt, non-asymptotic observers such that \v/v\ 
is constant perceive a thermal spectrum at temperature 



2tt 



(14) 



This is the case e.g. for (late-time) static observers at 
any radius r > r s : from the estimate Q, we compute 



= nv(r). 



(15) 



The time-derivative v(r) along a static trajectory can be 
read off from the metric (JTJ, as u(r) = (1 — r s /r) -1 / 2 , and 
therefore the temperature of outgoing modes measured by 
static UDW detector is readily found to be 



t (R)=T H (l-r s /r)- 1 / 2 . 



(16) 



Again, this is consistent with the Tolman equilibrium 
condition [5J. 



E. Non-static trajectories: quasi-temperature and 
adiabatic approximation 

It has been argued in [10] that the approximate con- 
stancy of T_(t) = \v-(t)/v-(t)\/(2tt), expressed as 



V- 



T_ 



I* 



«1, 



(17) 



is a sufficiently condition for the perception of thermal 
Hawking radiation. This "adiabaticity condition", con- 
sidered previously [2l[ in the context of evolving black 
holes, allows to study the response of UDW detectors 
along non-stationary trajectories [10J in a very straight- 
forward way, simply by computing T_ (r) . Heuristically, 
when ry_ <C 1, the response function 1Z(t,Q) of a UDW 
detector switched off (smoothly) at time r will be indis- 
tinguishable from a thermal spectrum with r-dependent 
temperature T1(t). In fact, even when r\ is not small, 
this relationship still holds in the ultraviolet limit, where 
£7>T_(t). 

Following this rationale, we will call quasi-temperature 
of a time-dependent spectrum 1Z(t, fl) a function T{t) 
such that, at any given time r, the detailed balance con- 
dition holds in the ultraviolet limit, 



11(t, -0) 



o n/T(r) 



K(t,SI) for |0| >T_(r). (18) 



See the Appendix for more details on the adiabatic ap- 
proximation. 



F. Two quasi-temperatures for Hawking radiation 

Our approach to the study of Hawking radiation along 
general trajectories is based on this adiabatic approxima- 
tion. For a given trajectory, we shall compute the quasi- 
temperatures of both outgoing and incoming modes as 



T±(t) = — 



v±(t) 



v±(t) 



(19) 



As we will see, these quasi-temperatures together with 
the corresponding adiabaticity parameters 



n- 



J2 



«1, 



(20) 



provide a valuable handle on the response of UDW de- 
tectors to Hawking radiation, especially in the large Vt 
limit. 



III. RESPONSE OF GEODESIC DETECTORS 

In this section, we apply the quasi-temperature for- 
malism to the case of geodesic trajectories in the 
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FIG. 4. The quasi-temperature of outgoing (continuous) and 
incoming (dashed) as a function of proper time along the inspi- 
ral geodesic (E,L) — (1,3) until it reaches the horizon. The 
inset shows the geodesic in the equatorial plane. The step 
at r ~ 300 corresponds to the moment when the trajectory 
enters the Hawking region. 



Schwarzschild region (v > 0) of the Vaidya collapse ge- 
ometry described above. 



A. Schwarzschild geodesies 

Timelike Schwarzschild geodesies are characterized by 
two orbital parameters: their energy per unit rest mass 
E and their angular momentum per unit rest mass L — 
r 2 ip. In Eddington-Finkelstein coordinates, the geodesic 
equations read 



r 2 + (l-r s /r)(l + L 2 /r 2 )-- 
-(1 - r s /r)v 2 + 2vr + L 2 /r 2 



E 2 



•1. 



(21) 



From these equations and (|3]), we can in principle ob- 
tain -u + = v, w_ and their derivatives — hence the quasi- 
temperatures T± and the adiabaticity parameters ij± — 
along any geodesic (E,L). This approach is illustrated 
in fig. H]for one inspiral trajectory with (E,L) — (1,3), 
by means of a numerical integration of (|21[) . This plot 
shows interesting phenomena: the setting off of Hawk- 
ing radiation after gravitational collapse, the increase of 
the quasi-temperature of outgoing modes, and of that of 
incoming modes as the horizon is approached. 

However, in order to get an analytical understanding 
of this plot, and in particular of the dependence of T± 
and r)± on E and L, the best approach is to consider spe- 
cial regimes of interest (circular orbits, radial trajectories 
in the asymptotic and near-horizon regimes) separately. 



This is the approach we shall take in the following. Our 
presentation will be rather pedestrian: we have found, 
indeed, that the details of the computation of T± and 
rj± provide a good deal of insight into the inner work- 
ings of the Hawking phenomenon; they display in a very 
transparent way the interplay between the geometry of 
spacetime (and its geodesies) and the structure of the in- 
vacuum which is responsible for the evaporation of black 
holes. 



B. Circular orbits 

Circular orbits are characterized by their Schwarzschild 
radius r > 3r s /2, or equivalently by 5 — r/r s — 1 > 1/20 
By virtue of their stationarity, these orbits have constant 
v, given by 



m 



5-1/2 



1/2 



(22) 



In particular, v = 0, hence T + = 0: incoming modes do 
not couple to UDW detectors on circular orbits. On the 
other, from (J3|) we get for the outgoing modes 



v-(S) 



kv(S) 1 + Se 



s _ KV W"(5e>-™) 



W'{5e s ' KV ) 



(23) 



Using (f2"2"j) and the standard formula for the derivative of 
the Lambert VT-function, 



W'{z) = 



W(z) 



we arrive at 



T^(v,5)=T H 



1 



z(l + W{z)Y 



1/2 



(24) 



1- 



8-1/2 
W(6e s ~ KV )(2 + W{5e s - KV )) 



[l + W{Se s - KV )Y 



(25) 



In this expression, the second term represents the tran- 
sient regime before the orbit has been "swallowed" in the 
Hawking region. After this transient, T_(t>, 5) reaches the 
stationary value 



Tl [l ' c {5)= lim T_(v,S)=T H 



1 



1/2 



1/2 



(26) 



1 The corresponding energies and angular momenta are given by 



E circ (r) = r s 



a 2 



(*+!)(* -1/2) 



1/2 



icirc(r) = r s 



(l + S) 2 



2(<5 - 1/2) 



1/2 



T-(5) 




1.0 1.5 2.0 2.5 3.0 



FIG. 5. Stationary trajectories: the temperature measured by 
static detectors (dashed line) and by geodesic detectors on a 
circular orbit (continuous line) around a Schwarzschild black 
hole as a function of S = r/r s — 1. 



Observe that, in spite of the fact that circular orbits 
have zero acceleration, this temperature is always larger 
than the temperature measured by static detectors at the 
same distance from the hole (compare with (|16p ). and di- 
verges on the photon sphere (r = 3r s /2). This can be 
interpreted as a transverse Doppler effect: as the photon 
sphere is approached, the relative velocity between cir- 
cular and static trajectories tends to the speed of light. 
These results are illustrated in fig. [5j 



C. Radial trajectories: surprises on the horizon 



Radially infalling trajectories (L = 0, r < 0, v > 0) are 
parametrized by their energy E, with E > 1 correspond- 
ing to unbound states and E < 1 to bound states. (The 
limiting case E = 1 describes a detector dropped from 
infinity with zero velocity into the black hole.) In addi- 
tion to the time-dependence of the Hawking phenomenon 
itself (the transient represented by w-dependent terms in 
(|25[) ). these trajectories are intrinsically non-stationary: 
as the black hole is approached (6 — > 0), the local Rie- 
mann curvature becomes larger and larger. To avoid mix- 
ing these two effects, from now on we shall assume that 
the trajectory is well inside the Hawking region, and focus 
on the time-dependent effects arising from the trajectory 
itself. 
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Outgoing modes. Consider first the outgoing modes. 
Using ([3]), we have 

(27) 



(5e s 



i-~ (Se s - KV )- 

_ 5(1 + S) + (5 - nv)(2S + 5(5 - kv)) - k5v 

5 + 5(5 — kv) 

Let us focus on the aforementioned limiting cases: the 

asymptotic limit (5 3> 1) and at horizon-crossing (5 = 0). 

In the asymptotic limit (which requires E > 1 to exist), 

the system (EJ) gives 5 (E 2 - l) 1/2 /r s , v ~ E- (E 2 - 

l) 1 / 2 , 5 ~ 0, and v ~ 0, hence from J37J with S > 1, 



K1> 



(/_ asymp 

Thus, in this limit we get 



-k(£+ \/£ 2 - 1 



T* symp (E) =T H [E+ VE 2 -1 



(28) 



(29) 



This formula is consistent with Hawking's original predic- 
tion, with a Doppler factor accounting for the relative mo- 
tion between the detector and the black hole when E > 1 . 
In this regime the adiabaticity parameter r/ a ^ ymp (E) is of 
course vanishingly small. 

At horizon-crossing (5 = 0), on the other hand, the 
equations of motion (|21|) give 5 = —E/r s and v = 1/2E, 
and, by differentiation of the radial equation, 5 = — l/2r 2 . 
Plugging this into ([2"7f now gives 



V- 

it- hor 

and therefore 



5 + 25(5 - kv) 



2E 



= -AkE 



T^ or (E) = 4ET H . 



(30) 



(31) 



Thus, not only is the quasi-temperature not zero on the 
horizon, but for unbound states (E > 1) it is actually 
larger than the Hawking temperature perceived by static 
observers at infinity. The formula (|31l) is consistent with 
the result of [10] 

We have stressed that the near-horizon regime is not 
stationary, hence that whether pip can be interpreted 
as a temperature depends on the value of r\- (E) there. 
Further differentiation of ([2~Tj) and (j2"7| gives 



/?_(£) 



1 
E 2 



(32) 



2 Considering an "observer freely falling from infinity", they find 
that "in the last stages of his approach to the horizon, and sur- 
prisingly at first sight, the effective temperature rises reaching 
exactly four times HawkingOs temperature" [HI , 



This number is of order 1 for all unbound trajectories 
(E > 1), and never smaller than ir/2 ~ 1.6: detectors 
with frequency ft ~ T^ or (E) will not confuse the vacuum 
state with a thermal state. Large-frequency detectors 
(Q > T hol (E)), on the other hand, will not be able to 
make this difference. 

Incoming modes. Let us now consider the quasi- 
temperature of incoming modes. As we saw earlier, these 
do not couple to UDW detectors on stationary trajecto- 
ries (static or circular); it may seem appropriate to as- 
sume that the same holds along radial trajectories. This 
is indeed the case in the asymptotic limit, where we 
have seen that v ~ — but what about the near-horizon 
regime? 

We have obtained v — 1/2E on the horizon. To get 
the corresponding value of v, we differentiate the second 
equation of (|2"Tj) . On the horizon, this gives v = —k/AE 2 , 
and therefore 



-rhor 



V ' 2E 



(33) 



For unbound states (E > 1), this quasi-temperature is 
always smaller than that of outgoing modes. The sit- 
uation is different for bound states (E < 1): for any 
E < 1/V8 we have 7* or (£) > T^^E), i.e. the incom- 
ing modes dominate the outgoing modes, and in the limit 
where E — > (detectors dropped with zero velocity from 
close to the horizon) , the quasi-temperature of incoming 
modes diverges. Taking these modes into account, we 
can therefore conclude that the slower it is approached, 
the "hotter" a black hole appears. In our view, this sur- 
prising result is the main contribution of this work. We 
will comment further on this finding in sec. IVII 

Again, the proper interpretation of this result should 
be derived from the value of the adiabaticity parameter. 
In this case, we get 



v\ OI (E) 



2tt\1-8E 2 



(34) 



Thus, just like for outgoing modes, r]^ 1 is bounded from 
below by a number of order one (actually 2ir ~ 6.2): that 
is, the adiabatic approximation is not good in this regime, 
and (|33l) can only be interpreted as the ultraviolet decay 
rate of detector spectra. 



D. Inspiral orbits: effect of angular momentum at 
horizon-crossing 

It is not difficult to generalize the above results to in- 
spiral trajectories with non-zero angular momentum. As 
before, the equations of motion ([21]) gives 5, 5, etc., as 
functions of E and L in the various regimes, and using 
it is a simpler matter to derive T^ >V (E : L). We do 



7\E) 



15 



10 





FIG. 6. Radially infalling trajectories: quasi-temperature of 
outgoing and incoming modes in different regimes, as func- 
tions of the energy E. For small energies E <C 1 (bound 
states), incoming modes dominate. 



not repeat this computation here; the results are 

T hor (E,L)=AET H (35) 



and 



rf"(£,L) = !| 



Sli- 



8£ 2 



1 + L 2 /r 2 s 



(36) 



and similar formulas for rf^ 1 which we do not give here. 
Thus, while the quasi-temperature of outgoing modes at 
horizon crossing is independent of the angular momen- 
tum of the trajectory, the quasi-temperature of incom- 
ing modes is not. Instead, it is given by the above non- 
monotonous function of E and L. 



IV. HAWKING FLUX ALONG RADIAL 
TRAJECTORIES 

It is well-known that the response of UDW detectors 
to vacuum radiation is a priori independent from the ex- 
istence of a non-zero energy-momentum tensor (T a &). In 
the Hawking effect (unlike the Unruh effect), however, 
it turns out that (T a b) ^ after the collapse: at future 
infinity, one finds (T uu ) = nTfj/12 (and all other compo- 
nents zero), confirming that a Schwarzschild black hole 
does produce a stationary outgoing flux — that it actually 
evaporates. 

What flux would infalling observers measure at 
horizon-crossing? Given the (surprising) results we have 
found for the response function of UDW detectors in the 



12F(E) 
nT 2 H 



-50 




FIG. 7. The flux J-(E) perceived by radially infalling ob- 
servers with energy E, both in the asymptotic region ("asym") 
and at horizon-crossing ("hor"). Compare with Fig. [6] 



limit E — > (also because they are based on an adiabatic 
expansion, which might seem uncompelling to some read- 
ers), it is interesting to consider the flux 



T(E) = (Ta„) u a n b 



(37) 



along radially infalling trajectories u a (E) as a function 
of their energy E, say in the direction n b orthogonal to 
u a in the outwards direction. This quantity has been 
considered previously by several authors, see |23 - t24| . Is 
this quantity also divergent as E — >• 0? 

To answer this question, we can rely on standard results 
on vacuum energy-momentum tensor in two-dimensional 
spacetimes: given null coordinates (v+,v-), the vacuum 
energy-momentum tensor can be written as [191 |25L [26J 



\J-v±v±) 



1Z7T 

RC 
( v ' v+) = 96^' 



where C is the conformal factor such that 



ds 2 = —C(v+,V-)dv+dv_ 



(38) 



(39) 



and R = — AC^ 1 d v+ d v _ InC the two-dimensional scalar 
curvature. 

In our case — which, by virtue of our assumption of 
spherical symmetry, is effectively two-dimensional the 
factor C (t>_|_, vJ) can be obtained from the Vaidya metric 
JTJ) by inverting the relation ([3]) between V- and (v,r). 
This gives 



C(v+,v-) 



-1 W( 



ae 



-a+Kv+ \ 



l + W(-ae- a+KV +) 



(40) 
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where a = 1 + kv_ 
considering the a —, 
we compute 



. Applying the formulas (|55|) and 
limit for kv 3> 1 (at late times), 



ric of the form 



\-^ V-V- 



\ hor 



_nTl 
12 

7rT 2 



(41) 



Note that the incoming component (T v+V+ ) at the horizon 
is minus the outgoing component (T v _ v _) at infinity. 

Contracting (|41l) with the radial 4-velocity u a (with 
components v = 1/2E and w_ = 2Ee~ KV at the horizon) 
and its outwards-pointing unit normal n b (with compo- 
nents n + = 1/2E and n_ = —2Ee~ KV ), we arrive at 



J* OT (E) = -ttT%(2E 2 



1 



48£ 2 



(42) 



This describes an ingoing flux, which — like the quasi- 
temperature T^ or — diverges when E — > 0. This expres- 
sion is to be compared with the flux at infinity, given 
by 



T asymp (E) 



nT; 



T li( E +^ 2 



1 



(43) 



We plot both functions of E in Fig. [7J Numerically, 
an observer dropped with zero velocity from e.g. r — 
100r s /99 (i.e. such that E = 0.1), the flux given by (g2J 
is ~ 25 times larger than the Hawking value ttT^/12. 



V. ROLE OF CURVATURE IN NEAR-HORIZON 
HAWKING RADIATION 

We have found that UDW detectors crossing the hori- 
zon of a Schwarzschild black hole will in general record 
non-zero quasi-temperatures as well as an ingoing flux. 
How much of this effect is due to the structure of the 
collapse vacuum, and how much to the mere fact that 
spacetime is curved on the horizon? It is impossible to an- 
swer this question within the framework of Schwarzschild 
black holes, because both their surface gravity (hence 
Th = k/2tt) and their Riemann curvature on the hori- 
zon are controlled by the same parameter, namely their 
mass. To get some insight into the role played by cur- 
vature in the Hawking effect, it is therefore necessary to 
consider artificial black holes with non-Schwarzschild ge- 
ometry, in the spirit of the "analogue gravity" program 



A. Artificial black holes with flat horizons 

Since we are interested in models of gravitational col- 
lapse, we shall consider generalizations of the Vaidya met- 



e(-v) + G{v)f(r))dv 2 + 2dvdr + r 2 dVt 2 (44) 



ds) 



where f(r) could in principle be essentially any func- 
tion of r. In particular, a function f(r) such that (i) 
f(r s ) = 0, (ii) f'(r s ) = 2k, (Hi) linv-,.^ f(r) = 1 
and (iv) linv^o f(r) = — oo will yield a black with the 
same asymptotic structure and surface gravity k as the 
Schwarzschild black hole, but with a different Riemann 
curvature distribution. Here we will focus on an artificial 
collapse metric of the form f(r) = (r/r s — l)^(r), where 
\&(r) ~ 1 near the horizon but ensures that conditions 
(m, iv) above are satisfied. This corresponds to a space- 
time where curvature (in the (v, r) sector) is concentrated 
on two locations: at the shell (v = 0), where it has a 6- 
function singularity, and at large radii (say r > r c , where 
r c is the extension of the flat region outside the horizon) , 
where the derivatives of ^>(r) start taking non-zero val- 
ues. So long as we only consider trajectories contained 
within the domain T = {v > 0, r < r c }, we can think 
of this model as describing a "flat black hole" . Does this 
spacetime radiate in the same way as the Swarzschild 
black hole? 

Repeating the computation leading to © in the Vaidya 
case, we find for this "flat black hole" 



.(»,r) = -2r s (l + Se- KV ). 



(45) 



Just like the Vaidya case, this equation defines a "Hawk- 
ing region" where v~ ~ — k" 1 ; the only difference dif- 
ference is this region grows relatively faster than in the 
Vaidya case, as can be confirmed by comparing the con- 
ditions 5e~ KV < 1 (flat) and 5e s ' KV < 1 (Vaidya). 



B. Quasi-temperatures at the flat horizon 

To estimate the quasi-temperatures measured by UDW 
detectors in geodesic motion around this flat black hole, 
we repeat the steps in sec. IIIII write the geodesic equation 
for trajectories with orbital parameters (E,L), which is 
now 



r 2 + /(r)(l + 
-f(r)v 2 + 2vr 



L 2 /r 2 ) = E 2 
+ L 2 /r 2 = -1, 



and from (J45|) compute 

5 — 2kvS — k5(v + kv 2 ) 



V- 



S — kSv 



(46) 



(47) 



Clearly, in the asymptotic limit, provided f(r) — > 1 for 
r — > oo, these equations lead to the same results as in 
the Schwarzschild case. This confirms that asymptotic 
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Hawking radiation is a global phenomenon which does 
not depend on the actual curvature distribution around 
the black hole. 

Consider however the case of UDW detectors in the 
near-horizon region, where f(r) ~ r/r s — 1 = 6. In this 
regime, (|46[) gives the same values for 6, v and S as in 
the Schwarzschild case, and in particular 5 = 2nv5, but 
now v = —kv 2 . Thus, all terms on the numerator of (I47p 
cancel, and therefore 



and therefore 



-hor 



(E,L) = 0. 



(48) 



The superscript indicates that this identity holds not 
just on the horizon, but in the whole near-horizon region 
where the black hole is flat. In this regime, the outgoing 
modes do not couple to UDW detectors. This shows that, 
for these to become thermal, there must be a region where 
the Riemann curvature is non-zero between the horizon 
and the observer (even if its actual distribution in space 
is irrelevant). 

As for the incoming modes, the relationship v = —kv 2 
gives 



Tl OI (E,L) 



T H (L 2 + r 2 ) 
2E r? 



(49) 



The incoming modes, which unlike the outgoing modes 
must have encountered some curvature on their way to 
the horizon, do couple to UDW detectors, even more so 
that L is large and E small. 

We have seen at the end of the previous section that 
there exists an "irreducible non-adiabaticity" of horizon- 
crossing radial geodesies in the Schwarzschild spacetime. 
This "irreducible non-adiabaticity" is also present in the 
fiat case, as can be confirmed by the computation of 77+: 
sticking for simplicity to L = 0, we get 



V ^(E) = 27T. 



(50) 



This irreducible non-adiabaticity is thus independent of 
the local geometry of the horizon; it is an intrinsic feature 
of Hawking radiation as perceived by freely-falling near- 
horizon observers. 



C. Flux at the horizon 

We close this section by repeating the computation 
of the flux measured by infalling observers J-(E) = 
(T a b)u a n b . In the fiat case, the relation (|45)) gives 

formula 



hence from the Davies-Fulling-Unruh 



F hoI (E) 



A8E 2 ' 



(52) 



Again, we hnd that flattening the horizon cancels the 
effect of outgoing modes at the horizon, but not that of 
ingoing modes. 



VI. DISCUSSION AND CONCLUSION 

In this paper, we have considered Hawking radiation 
from the perspective of UDW detectors evolving on non- 
asymptotic trajectories. When the trajectory is not sta- 
tionary, or just after the collapse, their response is not 
exactly thermal, but by using a "quasi-temperature" for- 
malism based on a suitable adiabatic expansion, we have 
shown how to get a basic understanding of their response, 
especially in the large frequency limit. This quasi- 
temperature formalism is especially convenient from the 
computational perspective: formula (III E|) is explicit and 
can be straightforwardly applied to any trajectory of in- 
terest. 

Following this approach, we have obtained several in- 
teresting results, which disprove the notion that Hawking 
radiation reduces to Unruh radiation at the horizon, or 
that "the equivalence principle is restored" [4] there (that 
is, assuming that the Hawking effect violates it in some 
sensejj). They are: 

• In spite of the fact that these trajectories are not 
accelerated, the temperature perceived on circular 
orbits is always higher than that on static trajecto- 
ries at the same distance from the hole — it actually 
diverges on the photon sphere. 

• In the near-horizon region, a freely-falling detec- 
tor couples to both the outgoing and the incom- 
ing modes of the field, to which it associates two 
different quasi-temperaturesQ depending on its en- 
ergy and angular momentum. Both can be arbi- 
trarily high (though in different regimes: for fast 
and slow moving detectors respectively) and nei- 
ther has a Dopplerian velocity-dependence relative 
to the hole. This singular behavior of Hawking ra- 
diation in the E — > limit is confirmed by a flux 
computation. 

• While the thermality of outgoing modes relative 
to near-horizon geodesic observers can be traced 



(T, 



V+V+ I 



12 : 



3 We probably would not have described the Hawking effect in 
those terms anyhow. 



(T v _ v _) = \T V _ V+ ) = (51) 4 At least in the s-wave approximation discussion in this paper. 
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back to the spacetime curvature, this is not the 
case for incoming modes: these appear thermal 
(both in terms of UDW response functions and of 
fluxes) to infallers also in spacetimes with flat hori- 
zons. In both cases, there exists an irreducible 
non-adiabaticity at the horizon, in the sense that 
|7±| ~ T\ at the horizon. 

In a nutshell: far from a no-particle vacuum, a detector 
dropped with zero velocity from near the horizon will per- 
ceive intense Hawking radiation coming from infinity — 
that is, from the sky. 



ACKNOWLEDGMENTS 

SS is supported by an SPM grant from the Council of 
Scientific and Industrial Research (CSIR), India. We are 
grateful to Thanu Padmanabhan for suggesting to study 
Hawking radiation from the perspective of infalling de- 
tectors and for kindly hosting one of us (MS) at IUCAA, 
India, during a research visit. We also thank Antonin 
Coutant, Ahmed Youssef and Carlo Rovelli for many use- 
ful conversations. 



Appendix A: Adiabatic expansion of UDW-like 
response functions 

The validity of the adiabatic approximation mentioned 
in sec. Ill El has been studied in some detail [13, l2a [29| 
in the context of non-uniformly accelerated UDW detec- 
tors in fiat spacetime (approximate Unruh effect). In this 
setup, Barbado and Visser [29| have devised an "adiabatic 
expansion" for the UDW transition rates in terms of the 
quantities a' n ' (r) / l a(r) n+1 , where a(r) is the instanta- 
neous acceleration of the detector. 

In this appendix, we describe the adiabatic expansion 
of any functional of a slowly varying function k(r) of the 
form 



R(T,rt;k] = / 

J — : 



dsK(ila)g(A T , s [k]). (Al) 



where K, g are any functions such that the above integral 
is well-defined, and A TiS [/c] is the functional 



A T , a [k] 



v(t) — v(t — s) 
v(t)s 



with v(t) any solution of the differential equation 
v(t) = —k[r)v{T). 



(A2) 



(A3) 



The UDW response functions of the kind discussed 
in this paper obviously fit this scheme, with K(z) — 



cxp(— iz), g(z) = ln(z), and k (divided by 2ir) the quasi- 
temperature function. 

The functional R(t, f2; k] defined by (jAlj) may be seen 
as a function of r, f2, and the infinitely many deriva- 
tives fc^(r) of k at a given instant r, or — better — as 
Q^ 1 times the dimensionless function 



R(v ,vx,...) 



duK(u)g(A T , u/n [k}) (A4) 



of the dimensionless ratios v n = k^(r) / T2 n+1 . If we 
denote i?o( J/ o) the value of R{vq,vi, . . .) in the partic- 
ular case where k(r) — k is constant, the question we 
wish to answer is: in the general case where k(r) is not 
constant, can R(yo,i>i, . . .) be systematically expanded 
about Rq^q), with each corrective term depending on 
only finitely many i/ n 's? 

The ke y id ea to get started with this question, which we 
found in [29( , consists in introducing a scaling parameter 

a and £;( q )(t — s) = k(r — as). Denoting z/„ = a n u n the 
corresponding scaled dimensionless ratios, we have 

1 d n 



R(fO,Vl,.. 



E; 

m=0 



i\da J 



-H[y , v x , 



) 



_ (A5) 

Using Faa di Bruno's formula (chain rule for higher 
derivatives), we may write the m-th term in this expan- 
sion as 



1 _ m . r°° 
-fJ2 duK(u)g^(A )B m JA 1 ^ 

m - j=0 J -°o 



• A m _j + i). 



(A6) 



where B n j are the Bell polynomials^, and 

Aj = diA TtU/n [k^]\ a ^ . (A7) 

In the next paragraph, we will show that Aj is given by 



A, 



^Sj(i/ u) 



(A8) 



for some function Sj. This result calls for several com- 
ments: 

• The m-th term in the Taylor expansion (|A5|) , given 
by (|A6[) . is is a function i? m (^o, • • • , v m ) of the first 
m + 1 parameters v n only, viz. 



R(i/ , v\,...) 



m=0 



Rm(vo, ■ • • , Vm). 



(A9) 



The Bell polynomials are defined by B m ^(z\,- 



2m-fc + l 



m! ^2 n^Li " t tt (^f") Ji i where the sum runs over all sequences 
of non-negative integers (ji, ■ ■ ■ ,.J m — fc+l) such that Ji + ■ • ■ + 
jm-h+1 = k and h + 2 J2 H h (m - k + l)j m - k +l = m. 
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• The relevant parameters controlling the conver- 
gence of the adiabatic expansion (|A9[) are indeed 
Vj/i>i +1 = fcW(r)/fe(rp' +1 , as observed in [H in 
the context of the Unruh effect; in particular, we 
have for each m > 1 



Ri 



.Oo,"- ,^m) =0( SUp -T^j-j 
^ l<7<m Vr\ ' 



(A10) 



When fc(r) is a polynomial function of r (so that 
only finitely many z/ n 's are non- vanishing) , we are 
assured that 

R{v ,vi,---) = Mvo)+o(sup-^). (All) 

In other words, the adiabatic approximation is good 
when all the expansion parameters Vj/vq are 
small. 

• The adiabatic expansion is not a high-frequency 
expansion; indeed, the functions 8j(vqu) = 
8j(k(r)u/£l) can well be non-polynomial in l/Q. 



To prove (|A8[) , we begin by expressing the argument of 
g in (|A4[) explicitly in terms of the dimensionless ratios 
/v n 's. From (|A3I) . we have 



K,u/n[k] = - da exp / dpfc(r + p). (A12) 

u J-u/a J a 

Expanding k{r + p) in powers of p and making the change 
of variable p — Qo~ in the er-integral, this gives 



" P n+1 »n 



Ar.u/nW - - / dji exp ^ f— i. (A13) 

w Vo fri (n + 1)! 



Next, we consider the effect of scaling by a on (|A17|1 . 

-| /"It °° 71+1 n 

u ^0 ^ (n+1)! 

and compute its j-th derivative with respect to a, as in 
(|A17|) . The 0-th derivative is easily found to be 



A 



e v ° u - 1 
v u 



and for, j > 1, 



-» - y^r^ +1 Ao. 



i + i 



(A15) 



(A16) 



Computing explicitly the derivatives of (|A15[) . we arrive 
at AMI) with 



(A17) 



(-1V+V ( ^ (-D r N 



z \ * — ' r\ 

7-=0 



We close this appendix by evaluating explicitly the first 
corrective term R\ in the particular case where K(z) = 
e~ lz and g(z) = lnz, namely 

/oo 
dse- ms ln(A TjS [fc]). 
-OO 

(A18) 
The real part of this integral is the UDW response func- 
tion along a trajectory with quasi-temperature T = k/2ir, 
studied in the main body of this paper. As is well-known, 
the 0-th order Rq(i>o) is the standard two-dimensional 
thermal spectrum 



-Ro(z'o) 



2tt 



e 2ir/v _ i ' 

From (fM|, dSJ and (|ATf|) , we have 

D i \ Vl 

Rx{vo,vi) = -oX 



(A19) 



°° due le^(l-^ + K U )V2)-l 

p"0« — 1 

-oo 



e "o" _ 1 



The integrand has simple poles at u n = ImnjvQ for each 
integer n^O. Closing the contour in the lower half-plane 
and summing over n < — 1, we obtain 

7/, / 7i" 2 (l + i7rcoth(7r/^o) 

^iK-i) = 4 V . , 2 , , : 7 | • ( A2I.) 

f \ <v smh (7r/^o) 

Reinstating the dimensionful variables fc(r), fc(r) and $7 
and writing T = k/2ir, this gives 



Re{i?(r,rj;fc = 27rT]} = 



2tt 



1 + (A22) 
(A23) 



n(e n / T - 1) 

T(r) g e a / 2T ( T ) 

T(r) 2 8ttT(t) sinh 2 (fi/2T(r)) " 



Remarkably, the second term in the brackets is exponen- 
tially vanishing in the ultraviolet limit, showing that the 
adiabatic approximation is excellent in this regime, what- 
ever the value of the parameter T(t)/T(t) 2 . 
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